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The derivation of the equation of one-dimensional movement of a solitary shock 
wave is given. This derivation shows, that the differential equation of movement 
of a solitary plane shock wave in the channel with variable area, is exact, if simpli- 
fying assumptions, made during derivation, are realized. But these assumptions 
in plane geometry it is possible to realize only approximately; situation with 
spherical and cylindrical shock waves is opposite. 



1 Introduction 



In 1957 Chisnell [T], being based on ideas of Chester [2], has deduced the 
equation of one- dimensional movement of a solitary shock wave on substance 
with polytropic equation of state. A little bit later the Whitham [3] has 
offered "... simple rule", which together with shock relations "... determines 
the motion of the shock wave". Actually Whitham was the first, who has 
considered the equation of one- dimensional movement of a solitary shock wave 
on substance with the arbitrary equation of state (see below the equation 
([14))) . but as the approximate equation, from which for a special case of the 
polytropic equation of a state, one can derive the Chisnell equation. 



2 The derivation of the equation 

We shall consider one-dimensional movements of solitary shock waves. Move- 
ment of a shock wave is identified with movement of its front. One-dimensional 
movements of shock waves are understood as movements of spherically sym- 
metric, spherically symmetric and plane shock waves. Shock waves are as- 
sumed solitary, i.e. it is not considered overtaking one shock wave by other 
shock wave etc. Naturally, the solitary shock wave is an abstraction. Be- 



low, during the derivation of the equation of one-dimensional movement of a 
solitary plane shock wave the differential calculus will be used, how it was in- 
tuitively used in Newton time. After development of the nonstandard analysis 
by Robinson (see jl]) such operating with infinitesimal differentials may be 
regarded as quite correct. 

Following Chisnell, we shall consider movement of a plane shock wave in 
the channel, having in some place infinitesimal jump of the sectional area dA 
- see figure (cp. it with Fig. 1 in Chisnell's article this figure can be also 
compared with Fig. 1, p. 193 in [£]). Let's choose the system of coordinates 




A A + dA X 



Figure 1: The shock wave, separating region 1, 2, is incident on small change 
in the area of a channel from A to A + dA. The resulting transmitted shock 
separates regions 4, 5. Regions 3, 4 are separated by simple Riemann wave. 
The shape of the channel is shown on the top of the figure. 

connected to the channel, which axis X we shall direct along the channel. We 
shall assume, that walls of the channel are indefinitely rigid, so it is possible 
to neglect interaction of a shock wave with them, the plane front of a shock 
wave is perpendicular to walls of the channel etc. (i.e. it is supposed, that 
it is possible to neglect the three-dimensional effects arising at movement of 
a shock wave), the substance before front of a shock wave is homogeneous 



and is motionless relatively the channel (at presence of the area change last 
assumption is implicitly containing in the assumption of homogeneity of sub- 
stance before front of a shock wave). As it is usually accepted, by an index "o" 
we shall mark parameters of substance before front of a shock wave. As one- 
dimensional movement is considered, all vectors can be identified with their 
projections to axis X. Thus: u = 0, p , po, £o, s - respectively, speed (= 
projection of speed on axis X), density, pressure, specific (referred to a mass 
unit) internal energy, specific (referred to a mass unit) entropy of substance 
before front of a shock wave. As the parameter, defining strength of a shock 
wave, we shall take: 

Po 

- the ratio of density of substance behind front of a shock wave to density of 
substance before front of a shock wave. It is universal, dimensionless, finite 
parameter. Chisnell as the parameter, describing strength of a shock wave, 
have used the ratio of pressure: z = p/p , - that can not be used, for example, 
for some equations of a state of substance, when p = 0. In one-dimensional 
case assignment of a (or z), with known parameters of substance before the 
front of the shock wave uo = 0, po, Po completely determines parameters of 
substance behind front of a shock wave and speed of movement of shock wave 
D - see theorem 5.5 in j^J. 

On Fig. 1 the situation is schematically represented, when the shock wave, 
moving from the left to the right in positive direction of axis X (D > 0), will 
pass through jump of area of the channel dA. Arising gas-dynamic break splits 
into simple Riemann r-wave (according to definition in § 16, [5]) moving back 
from the front of the shock wave and shock wave with changed on da strength, 
moving to the right (see § 17 in |3]). For initial shock wave moving from right 
to left (D < 0) instead of r-wave we should have Riemann /-wave - theorem 
16.2 in j^J. Situation, represented on Fig. 1 corresponds to supersonic motion 
behind front of the shock wave \u\ > c, where c - speed of a sound behind 
front of a shock wave. Following reasonings are correct without change and in 
a subsonic case - \u\ < c. 

Changes of speed u and pressure p behind front of the shock wave, corre- 
sponding to change of shock wave strength after passage of the area change, are 
" compensated" by changes of speed o^w and pressure dxp at adiabatic flowing 
of substance behind front of shock wave through the jump of the sectional area 
and by changes of speed and pressure in Riemann wave (for Riemann r-wave, 
when D > 0, index " + " will be used, for /-wave - index "_"). Thus, we have 



the system of the equations: 



du 

— ■ da = dAU + d±u, (2) 
da 

■ da = dAP + d±p. (3) 

da 

At adiabatic flowing of substance behind front of the shock wave through the 
jump of the sectional area dA following relations must be satisfied: 

d A (Apu) = 0, (4) 

^(t +£+ p) =0 ' (5) 

d A s = 0. (6) 

First of them expresses law of conservation of mass, the second is Bernoulli 
equation, the last - adiabatic condition. The adiabatic condition implies equal- 
ities: 

d A p = c 2 ■ d A p, (7) 
P 

d A e = — ■ d A p- (8) 
p 2 

Using (J7J) and (jHJ), from (jU)-© we receive system of the equations: 

d A u = -— ■ (9) 
u p 

d A p = c 2 ■ d A p, (10) 

The system of the equations (|9*|)- (fTT|) completely determines changes of param- 
eters of substance at its motion through the jump of the sectional area of the 
channel. 

In Riemann wave the changes of parameters of substance are in accord 
with relations: 

d ± u±^ = 0, (12) 
pc 

d±p = c 2 ■ d±p. (13) 

The equation (J12)) follows from definition of Riemann invariants and the equa- 
tion (|13j) follows from adiabatic condition for movement in Riemann wave. 



Using ©-dJ and (US)-®, from ©-© we receive the equation: 



.^ = (l ± l).(p. ± L.^). dtr . (14) 

A \u c J \da pc da J 

For a shock wave, moving from left to right (D > 0), in the right part (jTj| it 
is necessary to take signs " +" . If D < (the shock wave moves from right to 
left), in the right part of the equation (JT4"j) it is necessary to take signs "— ". 
It is interesting to note, that the right part of the equation (jl4j) is a product 
of a combination of speed of substance and speed of a sound on differential of 
appropriate Riemann invariant. 

For the channel with arbitrary (smooth) dependence of the area of section 
on coordinate x - A(x), in the assumption of plane symmetry of considered 
movement, it is possible to pass from differential (fT4"j) to integrated equality: 

A a 

dA' f fl 1\ f du 1 dp\ , , 

11 ±T -hj3±--^ - da', (15) 



Ai 



A' J \u c J \da' pc da' 



- where Ai and Oi - values A and a at some moment of time tj. 

The equation (JT5|l can be named (integrated) equation of movement of a 
solitary plane shock wave on homogeneous substance in the channel of variable 
area, as at known dependence A (x) it determines speed of shock wave D as 
implicit function of coordinate x and with trivial equation: 

f = jD(x) ' (16) 

- allows to describe movement of a shock wave completely. The equation 
(|14|) now, accordingly, can be named the differential equation of movement 
of a plane solitary shock wave on homogeneous substance in the channel of 
variable section. 

The equation ()14j) can be trivial generalized, if substance before front of 
a shock wave is not homogeneous. For this purpose, in equations (J2J), © it 
is necessary to replace differentials ^ • da and -# • da with the more general 
expressions du and dp, because using theorem 5.5 in [5] one can write: 

du du du du 

du = — -da+ — ■ du Q + — ■ dp Q + — ■ dp , (17) 
da duo dpo dpo 



(18) 



- instead of a, u , p , p , other choice of parameters is possible also. After 
that, from (|4"j)- (jl3|) . we can finally receive following differential equation of 
movement of a solitary shock wave: 



It is supposed, that appearing in (fTTjl - fTDJ) functions A(x), u (x), po(x), 
Po (x) are smooth. For the first time similar generalization was proposed 
by Vakhrameev - see [S] , [7] . 



3 Chisnell's solution 



Chisnell in article pQ has derived analogue of the equation (J 14)) for a special 
case of polytropic equation of state: 



j_ _ dA{z) 
A dz 



~fZ T (z-1) 



(7+1) 



2{( 7 +l)z+( 7 -l)} 



+ 

X 



72:{(7-l) Z +( 7 +l)} 

1 _ ( 7 +l)-(2-l) 



X 



(20) 



2{( 7 +l)z+(7-l)} 



(7-l)z+(7+l) 
2(z-l) 



Here 7 is adiabatic exponent for polytropic gas, z = p/po- The parameter z, 
describing force of a shock wave, is determined correctly, as for the polytropic 
equations of a state po 7^ implies po 7^ 0, if only the temperature is not equal 
to absolute zero. 

The indefinite integral of the equation ()20j) found by Chisnell, looks like: 



Af(z) 



constant. 



(21) 



Where 



/(*) 



Z7 ( z -l).l z + £L 



2( 7 -l) 



Let's look, how in case of polytropic equation of a state the equation (j2*U|) can 
be derived from the equation ([14)1 . For definiteness we shall consider shock 
wave, moving in a positive direction of axis X (D > 0). Then we take signs 
"+" in (HU): 



1 dA 

A dz 



1 1 
- + - 

u c 



f du 1 dp 

\ dz pc dz 



(23) 



Parameters of substance behind front of a shock wave are determined by 
Rankine-Hugoniot shock relations, which for polytropic equations of a state 
can be written in the following form: 



u{z) = {z-l) 



P (*) = Po 



2p 



c(z) = 

Differentiation of ()24)1 gives: 

2p 




Po {(7 + 1) z + (7-l)} 

(7 + l)z + ( 7 -l) 
( 7 -l)z+( 7 + l)' 

7^Po _ (7 ~ 1) g + (7 + jj 
Po ' (7 + l)-z + (7-1)' 



po {(7 + 1)2 + (7-1)}. 



2{( 7 + l)z + (7-l)} 



(24) 
(25) 
(26) 



(27) 



Substituting p4]l - (j?7jl in flZfl), we receive Chisnell equation (jUJ). 



4 Comparison of different solutions 

One of the most important examples of one-dimensional movements of shock 
waves is movement on homogeneous substance solitary converging spherically 
symmetric shock wave. Choosing small enough solid angle with top in the 
centre of symmetry of a shock wave, we can see, that it is possible to consider 
movement of a shock wave inside such solid angle, as movement in the channel 
with the variable sectional area (^4 = constant ■ r a , a = 2, r - radius of 
front of a shock wave in spherical system of coordinates with the beginning 
in the centre of symmetry of a shock wave), and all assumptions of section 2 
are precisely carried out. Therefore movement of spherically symmetric shock 
wave in spherical system of coordinates, which beginning coincides with the 



centre of symmetry of a shock wave, is described by the equation (see the 
equation (fT5j) ): 



— a ■ In— 



I 



( 



1 1 
- ± - 

u c 



) 



( 



du 1 dp 
da' pc da' 




(28) 



Similar reasonings show, that the equation ()28|) describes also movement of 
solitary converging cylindrically symmetric shock waves (a = 1, r - radius of 
front of a shock wave in cylindrical system of coordinates with the beginning 
on an axis of symmetry of a shock wave). Signs in the right part of the equation 
fl2~%|) are taken according to agreements in Sec. |2J). 

The finding of the analytical solution of the equation (J28|) for the concrete 
equation of state of substance can be not trivial problem - see Chisnell's 
solution (}2~T|) - (j22|) . However, with the help of asymptotic analysis of functions 
of real variable (see, for example, chapter 5 in 8J) it is often easily to receive 
asymptotic solution (at r — > 0) of the equation (|2~%j). that is also interesting 
enough. Thus received asymptotic solution of the equation (|2*%j). generally 
speaking, are logarithmically equivalent to the true solution - see definition 5, 
§ 1, chapters 5 in 8J. It is dictated by structure of the right part of the equation 
f!28|) . Analytical Chisnell solution for polytropic equations of state allows to 
receive strongly equivalent solution of the equation (|2*Hj) (see definition 4, § 1, 
chapters 5 in jH]): 



Here pi - pressure behind front of the shock wave, which is taking place at 
some moment of time U on radius (it is supposed, that pi > p ); a - is 
determined above, v is expressed by the formula: 



In the table exponent av are given for 7 = 5/3, 7/5, 6/5 together with the 
corresponding values from self-similar solutions: Guderley, Butler - see [T], P, 
[TU] . Landau and Stanukovich - see § 64 in jTTj, see also [T^ . 

It is visible, that the accordance is very good. Small distinction of values 
could be tried to explain, how it tried to make Chisnell P, that at derivation 
of the equation (fTlj) possible change d 2 a of strength of shock wave, caused 
by reflection of the simple Riemann wave, which has arisen at jump of the 
sectional area of the channel dA 2 , from earlier arisen distortion (at jump of 




(29) 




(30) 



Table 1: Values of the module of exponent according to calculations of different 
authors. 





Cylindrical wave (a = 1) 


Chisnell 


Butler 


Guderley 


Landau, 
Stanukovich 


7=1 


0,326223 


0,322441 






7 = 1 


0,394141 


0,394589 


0,396 


0,398 


7 = 1 


0,450850 


0,452108 








Spherical wave (a = 2) 


Chisnell 


Butler 


Guderley 


Landau, 
Stanukovich 


7 = 1 


0,652447 


0,641513 






7 = 1 


0,788283 


0,788728 




0,789 


7 = 1 


0,901699 


0,905385 







the sectional area of the channel dA\) behind front of a shock wave. But d 2 a 
must have second infinitesimal order, because d 2 a is bilinear function of dAi 
and dA 2 - d 2 a « dA\ ■ dA 2 (d 2 a can have first infinitesimal order for divergent 
shock waves). 

To understand the true reason of discrepancy of solutions we shall return to 
derivation of the equation (JTJJ). The derivation of the equation (fTlj) is unusual. 
It is not local in the sense, that changes of the values relating to various points 
of space are considered. Therefore it would be possible to expect, that the 
type of symmetry of task will somehow show itself. Equation (|T2*j) determining 
linear connection of infinitesimal changes of speed of substance and pressure in 
simple Riemann wave in case of spherical (cylindrical) symmetry is incorrect. 
It would be necessary to replace it with more general relation (see, for example, 
H3|, § 2.7 of chapter 2): 



1 acu 
D±u ± — D±p ± = 0. 

pc r 



(31) 



Where 

D - = m + iu±c) -^ (32 > 

a = 0, 1,2 - correspondingly, in case of plane, cylindrical and spherical sym- 
metry Thus, the equation (fTIj) for spherical (cylindrical) sound waves should 
be considered as approximated. 



5 Conclusions 



Given derivation of equation of movement of a solitary shock wave in the 
channel of variable area shows, why this equation may be only approximate 
for cylindrical and spherical convergent waves. But such approximation has 
high precision and may be used for estimations. 
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